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Abstract 

^ \ 
O . 

\ Two phenomenological models describing an SU{N) quark-gluon plasma are 

J ^ I presented. The first is obtained fi'om high temperature expansions of the free 

^ • energy of a massive gluon, while the second is derived by demanding color 

■ neutrality over a certain length scale. Each model has a single free parameter, 
, exhibits behavior similar to lattice simulations over the range Td — 5Td, and 
[ has the correct blackbody behavior for large temperatures. The N = 2 decon- 

CN ' finement transition is second order in both models, while N = 3,4, and 5 are 

' first order. Both models appear to have a smooth large- limit. For > 4, 

■ it is shown that the trace of the Polyakov loop is insufficient to characterize 
^ . the phase structure; the free energy is best described using the eigenvalues of 

I the Polyakov loop. In both models, the confined phase is characterized by a 

I mutual repulsion of Polyakov loop eigenvalues that makes the Polyakov loop 

' expectation value zero. In the deconfined phase, the rotation of the eigen- 

Q^' values in the complex plane towards 1 is responsible for the approach to the 

^ . blackbody limit over the range — bT^. The addition of massless quarks 

D i in SU{3) breaks Z(3) symmetry weakly and eliminates the deconfining phase 
transition. In contrast, a first-order phase transition persists with sufficiently 
heavy quarks. 
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I. INTRODUCTION 



The equation of state for the quark-gluon plasma is of great interest in several different 
areas of physics. In heavy ion physics, astrophysics, and cosmology, the equation of state 
is needed as input. There are two first-principles methods to obtain the equation of state: 
perturbation theory and lattice simulations. Only lattice gauge theory techniques are able 
to determine the equation of state directly over all temperatures but researchers in other 
fields have typically contented themselves with the extraction of a few important parameters 
from lattice results. When the equation of state is needed, the bag model equation of state 
|Jl[] is very often used, even though it is a poor representation of lattice results. 

As better alternatives, we have developed two simple models for the SU (N) gluon plasma 
equation of state. Both models are obtained by combining simple phenomenological ideas 
with Z{N) symmetry and the known features of the perturbative equation of state. These 
models have many desirable properties, and exhibit thermodynamic behavior similar to that 
obtained from lattice simulations. In particular, they give a reasonable picture of the crucial 
region from the deconfinement temperature to 5Td which has not been obtained by other 
means. The phase transition in these models is second order for N = 2, and first order for 
iV = 3, 4, and 5. It appears that the large- limit will be smooth in both models. Quark 
effects are easily included. 

In both models, the eigenvalues of the Polyakov loop determine all thermodynamic prop- 
erties. The Polyakov loop is defined for Euclidean finite temperature gauge theories as 



where Tis the temperature, and T denotes Euclidean time ordering. The Polyakov loop 
is is the natural order parameter of the deconfinement transition in pure gauge theories. 
Its trace in the fundamental representation can be related to the free energy F of a heavy 
source: 



where the expectation value is taken over a thermal ensemble of states. It has been known 
for some time that {Trp P{lt)) is of fundamental importance in describing the deconfining 
phase transition p| . In a pure gauge theory, the action is invariant under a global symmetry 
associated with the center C of the gauge group, but the Polyakov loop is not invariant. If z e 
C, then the symmetry requires z {Trp P{'^)) = {Trp P{'^)) , which implies {TrpP{lt)) = 
0. This is in turn interpreted as F = oo, indicating confinement. In the phenomenological 
models developed here, the eigenvalues of the Polyakov loop are the essential degrees of 
freedom rather than {Trp P{lit)) alone, a possibility also recently explored by Pisarski [Q. 
The Polyakov loop P is unitary for SU (N) gauge theories; we will denote its eigenvalues in 
the fundmental representation by 



after a diagonalizing unitary transformation. The eigenvalues in other representations will 
be linear combinations of the phase factors exp{i6j). For simplicity, we will refer to both 
exp(i6'j) and 6j as eigenvalues, relying on context to differentiate them. 




(1) 



^-F/T ^ l^TrpP{lt)) 
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Pjk = exp {i9j ) 5jk 



(3) 
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In our models, we make a mean-field assumption that the Polyakov loop eigenvalues 
are constant throughout space, and assume that the free energy / is a function of the 
eigenvalues. Confinement is obtained from a set of eigenvalues which make {Trp P{lt)) = 
0. As we show below, this is naturally obtained by a uniform distribution of eigenvalues 
around the unit circle, constrained by the unitary of P. In other words, confinement at low 
temperatures is a consequence of eigenvalue repulsion. In a pure gauge theory below the 
deconfinement temperature T^, the eigenvalues are frozen in this uniform distribution. As 
T moves upward from T^, the eigenvalues of the Polyakov loop rotate towards ^ = or one 
of its Z{N) equivalents, and P moves towards an element of the center. In the case of a 
first order transition, the eigenvalues jump at T^. In the models developed below, it is this 
motion of the eigenvalues which is responsible for the approach to the blackbody limit over 
the range — ST^. 

The next section briefly derives the conventional one-loop expression for the free energy 
of gluons in a constant Polyakov loop background. Using this expression, section |ITI| explains 
in detail why a naive Landau-Ginsberg treatment of deconfinement based on Trp P{lt) as 
a single complex order parameter fails for SU{A) and higher. Sections ^ and [V| introduce 
model A and model B, respectively. Model A is particularly tractable analytically in the 
cases of SU{2) and SU{?>). Section |VT| presents results for the pressure p, the internal 
energy £, and the interaction measure A for both models for = 2, 3, 4, and 5. Section |VI1| 
considers the effects of massive and light quarks in model A for 5'f/(3). A final section 
discusses our results. 



II. PERTURBATIVE EOS AT HIGH T 



There are several reasons for beginning with the perturbative expression for the free 
energy. As a consequence of asymptotic freedom, the perturbative expression for the free 
energy as a function of the Polyakov loop eigenvalues will be valid at sufficiently high tem- 
peratures. Any purely perturbative calculation will give a free energy of the form 



U,rt = T^F{P,g{T)) 



(4) 



where F depends directly only on dimensionless variables: the Polyakov loop P and the 
running coupling constant g{T). The leading order behavior, which is sufficient for our needs, 
is independent of g{T) and gives the blackbody behavior expected at high temperature. In 
addition, the perturbative expression for the free energy will give us significant insights into 
the structure and interpretation of possible terms in the free energy. 

At one loop, the free energy for gluons in a constant Aq background can be written as 



UUd) = In [det (-D 



)2 

adj 



where Dadj is the covariant derivative acting on fields in the adjoint representation, 
can be written in an unregularized form as 
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where the trace is in the adjoint representation and /3 = 1/T. Using a standard identity 0, 
this can be written in a famihar form as the sum of a zero-temperature contribution, which 
is divergent, and a finite temperature contribution: 

W) 



j^^trA ^-Uk + In [1 - exp {-{Suu - ^/3v4o)]| 



(7) 



where uik= k and the factor of 2 comes from the sum over polarization states. Disregard- 
ing the zero temperature contribution, the free energy density can be written as 



N 



fpert{e) = -Y: 2(1 



j,k=l 
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where the factor involving the Kronecker 5 projects out the singlet state. We have defined 
the differences of the fundamental representation eigenvalues as /\6jk = Oj — Ok] they are 
just the A^^ — 1 eigenvalues in the adjoint representation. It is instructive to expand the 
logarithm: 

N 1 . ^3; 



1 ^ r d-^h °° 1 

/ = -- y 2(1 - / y - 



(9) 



Each factor of n is associated with paths that wrap around space-time n times in the Eu- 
clidean time direction 0. The integral over k can be performed, yielding 

o N 

j,k=i 



fpert (0) — — 



N 



JnAe 



'jk 



71=1 



(10) 



Note that /„ is real because A6jk = —A6kj. The summation over n can be done exactly. 



giving 



fpert (0) 
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where |A6'|2^ = |A6'| mod 271 axid B^[x) 



This leads immediately to 



N 
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UUd) = y (1 - j-S,k) 

j,k=l 
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2x^+x'^ — -^ is the fourth Bernoulli polynomial. 



487r2 



jkh-w 



- 27r 



(12) 



which reduces to the usual black body formula in the case Aq = 0. The free energy is a 
smooth polynomial in the A6 variables, except at the points A6 = 0, ±27r, ±47r, ... where 
the mod function acts to maintain the periodicity which is manifest in the original form of 



The minimum of fpert occurs at = and values related to = by Z{N) symmetries 
for all values of the temperature T. We can characterize these solutions as 

e, = ^ n = 0,..,(iV-l). (13) 

The one-loop expression for the free energy of gluons propagating in the background of a 
non-trivial, but constant, Polyakov loop thus predicts a gas of gluons would always be in the 
deconfined phase. There is no indication that higher orders in perturbation theory modify 
this result 0. 
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III. FAILINGS OF THE LANDAU APPROACH 



We define Pp to be the expectation value of tlie trace of tlie Polyakov loop in the funda- 
mental representation, (Trp P{lt)) . Because Pp is the order parameter for the deconfining 
phase transition in pure gauge theories, it is natural to attempt to construct a Landau theory 
for the deconfining phase transition as a polynomial in Pp [p^p . 

An obvious form for the Landau free energy fp is 

fp = a,P*Pp + I {P^Ppf + I {P^Ppf + ^-^{P^ + Pf) (14) 

where the coefficients are all real and temperature-dependent. Most of the terms in fp 
are invariant under a U{1) symmetry Pp —>■ e^°'Pp. The term b]sf{PP + Pp^) breaks this 
symmetry down to Z{N), the correct symmetry group for the SU{N) deconfining transition. 
An exception occurs for SU (2): Pp is real, so the 62 term can be absorbed into 02, and fp has 
only a Z{2) invariance. While additional terms can be added to fi, this is the minimal form 
necessary to reproduce all known behavior in 3 + 1 dimensions. Models of the lattice SU (N) 
deconfinement transition valid in the strong-coupling limit indicate a first-order transition 
for all iV > 3 [0,0]; there is also an argument based on Schwinger-Dyson equations that 



the N = 00 limit also has a first-order transition |jTj]. Recent lattice simulations of S'f/(4) 
at finite temperature provide very good evidence that the SU{4:) transitition is first order 
p!5| . The deconfining SU{N) gauge transition can be first order for all > 3 with the 
simple Landau free energy fp only if the coefficients are carefully chosen functions of the 
temperature. In the large A^ limit, these coefficients must scale in such a way that fp oc A^^. 

There are several problems associated with using to obtain an equation of state over 
a wide range of temperatures. First, there are generally at least four unknown coefficients 
to be determined. Because Pp is dimensionless, the coefficients are all of dimension four in 
3-1-1 dimensions. We can safely assume that they grow no quicker than for large T, 
in accord with the usual black body behavior. The standard Landau assumption would be 
that each coefficient is a polynomial in T of degree four or less. This gives us generically 20 
undetermined parameters. Without additional assumptions, it is impossible to make useful 
progress. 

Secondly, this form of the free energy includes only part of the known high-temperature 
physics. As we have seen in the previous section, the high-temperature partition function 
depends on terms of the form Tta [P"') where n can be associated with a gluon trajectory 
winding around space-time n times in the temporal direction. Because Trp (P) Trp (P)* = 
TrA (P) + 1, the 02 term is associated with the n = 1 gluon trajectories, and so on. We see 
that a polynomial in Pp does not take into account all of the high-temperature perturbative 
result. 

Finally, the most general form of the free energy does not depend on P solely through 
Pp and Pp. The high-temperature perturbative free energy fpert illustrates this point. It 
is not a function solely of Pp and its conjugate, and cannot be written as an infinite series 
in Pp and Pp. A simple example will illustrate this point. Consider two diagonal matrices 
lying in S'f/(4), defined by 
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i7r/2 



-i7r/2 



(15) 



e-m/2 y 



and 



t/2 



V 



3i37r/4 



3i57r/4 



(16) 



i77r/4 ^ 



Both C/i and U2 have zero trace in the fundamental representation, yet the traces of their 
square are differenct: TrpUf = —4 and TrpU2 = 0, estabhshing that TrpU'^ cannot be a 
function solely of Trp U. Explicit calculation shows furthermore that /pert (Ui) 7^ /pert (^2)- 
At first sight, this may seem to contradict two standard results: a) the characters form 
a complete, in fact orthogonal, basis on class functions; b) all characters may be obtained 
from the fundamental representation by repeated multiplication and the application of 



Xa{U)xb{U)=Y:n{a,b;c)xc{U) 



(17) 



where all n's are non-negative integers. Taken together, these results might suggest that all 
group characters are polynomials in TrpU and its complex conjugate. Consider, however, 
the product representation N ®N . It is reducible into A'"(A'" + l)/2 © N{N — l)/2, which are 
symmetric and antisymmetric representations, respectively. Their characters are respectively 



Xs 2 
1 
2 



Xa 



{TrpUf + TrF (u 
{TrpUf - Ttf iy^) 



(18) 



Note that the sum xs + Xa is a polynomial in TrpU, but xs and xa are in general not. For 
example, in SU{3) the restriction on the eigenvalues of U imposed by det {U) = 1 allows us 
to prove 



1 r 



{TrpUf - Trp (c/^)] = Trp (u^) 



(19) 



in accord with the SUiZ) result 3 3 = 6 © 3. In this case, it is true that Trp (U^) can 
be written as a polynomial in TrpU and TrpU^ ■ However, for SU{A), unitarity of U gives 
instead 



{TrpUf - Trp [u^)] = - [{TrpUf - Trp (u^) 



(20) 



which shows that the 6 representation of SU{4:) is real. SU{N) group characters can be 
written as polynomials in Trp (U) and its complex conjugate only for SU{2) and SU{3). 

An alternative statement is that the Polyakov loop in the fundamental representation, 
TrpU, is not sufficient to determine the eigenvalues of U, beginning with the case of SU{4). 
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Let us label the eigenvalues for the Polyakov loop as Zi, with i = 1 to N. For SU{3), the 
characteristic polynomial for the Polyakov loop is 



i=l 



Z^ - {Zi + Z2 + Z^i) Z^ + {ZxZ2 + Z2Z-i + 2:3^1) Z - Z^Z^Z-^. 



(21) 



Since the determinant of a special unitary matrix is 1, we have -21-22-23 = 1, and the charac- 
teristic polynomial is 



z^TrpU + zTrpU- 



(22) 



so for SU{3), knowledge of TrpU determines all eigenvalues, and the free energy can be writ- 
ten as a function of TrpU alone. For SU{A), similar considerations allow the characteristic 
polynomial to be written as 



1 



/ - z^TrpU + z^- [{TrpUf - Trp (u^)] - zTrpU^ + 1 



(23) 



and knowledge of TrpU must be supplemented by the value of Trp (f/^). As increases, 
more information must be supplied to reconstruct the eigenvalues. 

IV. MODEL A 



We want to add terms to fpert which are subleading in comparison to the behavior of 
fpert, and introduce a mass scale into /. This mass scale will determine the deconfinement 
temperature T^. The first of our two models is obtained by adding, by hand, a mass to 
the gauge bosons, and working with the high temperature expansion of the resultant free 
energy. 

As before, we have 



/ 



1 ^ 



A^ 



(2^ 



E 

n=l 



-nl3wf,+in^9 jf. 



n 



(24) 



but now ujk = vA?-|-M^. It is easy to derive the first two terms in the high-temperature 
expansion. Higher order effects, which include terms of order TM^ and In (T/M), can 
be derived by more sophisticated methods. Such terms were derived in the case of a trivial 
Polyakov loop by Dolan and Jackiw |jl6[. See 0] for a relatively simple derivation of the 
general case. The result is 



fA{e) 
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E 

j,k=l 



TT 



'AO 



jk 



2^2 



2ti 



-Bo 



'AO 



jk 



2n 



(25) 



where the second Bernoulli polynomial B2{x) is given by B2{x) = x^ — x + 1/6 on the 
interval (0, 1). The above expression defines the free energy for our model A. In fact, the 
full one-loop free energy for a massive gluon always favors Aq = 0. We stress that we view 
this derivation as merely providing an indication of the type of additional terms that might 
appear in the complete free energy. The Bernoulli polynomials appear naturally in class 
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functions which are ahiiost everywhere polynomials in the ^'s; these class functions are thus 
well-suited for the construction of a Landau theory in the eigenvalues. In more explicit form, 
/a is given by 



N 



fA (9) = -T' E (1 - 



j,k=i 



N 



1 



+TV ^ (1 - -6,,) 

j,k=l 



(26) 



Because the free energy density Ja is a class function of the Polyakov loop eigenvalues by 
construction, Ja is invariant under gauge transformations. This important property will also 
hold for model B. 

The T'^ term has a clear minimum when A6jk = 0, corresponding to = and its 
Z{N) equivalents, and will dominate for large T. The term, however, will dominate 
for T small and has a global maximum at AOj/,. = 0. To explore the phase structure, 
it suffices to consider real and positive Pp; it is convenient to use a parametrization for 
the SU{N) angles which we will refer to as the parametrization. For N even, we rep- 
resent diagonal matrices as diag 
ordered such that tt > 4'n/2 > •• 
diag 
that 



exp (i(f)N/2, ■■, i4'ii ~^0Af/2 j ! with the eigenvalues 

> 01 > 0. For N odd, we represent diagonal matrices as 

exp {i(f)(N-i)/2, ••) 0, — .., —i(t>{N-i)/2j , with the eigenvalues again ordered such 

> <t>{N-l)/2 > •• > 01 > 0. 



TT 



At low temperatures, the 0will minimize 



N j 
j=2 k=l 



(27) 



For even N, this reduces to 



N/2 



^[4iV02_87r (2j-l)0, 



(28) 



The minimum occurs at 



(f)j = 



7r(2j-l) 

N 



(29) 



which is precisely uniform spacing around the unit circle. For odd N, we have a similar 
reduction to 



{N-l)/2 

^ [4iV02 - I67rj0, 



(30) 



and the minimum is given by 



(f)j = 



27rj 

N ' 



(31) 
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Note that in the case of even, unitarity forces the eigenvalues away from 6 = 0, e.g., for 
= 4, the four angles form an "X" rather than a "+". 

In this model, the deconfinement transition arises because of competition between the 
T'^ term, which tends to force all eigenvalues to zero, and the term, which forces the 
eigenvalues apart. We next discuss the analytically tractable cases of SU{2) and SU{3). For 
SU{4) and higher, model A is conveniently solved by numerical methods. 



A. Model A for SU{2) 

In the case of SU (2), we have only 0i = and 0_i = —0. At low temperatures, we have 
= 7r/2, giving Pp = 0. The free energy is 



2rp4 4^4 ^2rp2 ^2rp2 

^ = -^ + ^0^(0-vrf + ^ + ^0(0-7r) (32) 



This equation has an obvious symmetry under (f) ^ tt — (p associated with Z{2) invariance. 
It is convenient to define a new variable ip = n/2 — (p, which better manifests the Z{2) 
symmetry. We obtain 

/ = + -(TT^-Aip^) + — (7r2-4^2) (33) 

15 127r2 V J 4 47^2 V ^ J ^ ^ 

with ip = representing confinement. The phase transition is second order, in accord with 
the universality argument of Svetitsky and Yaffe f^. The deconfinement temperature is 
given by Ta = (3/2)^^^ M/tt ^ 0.38985M, and ^ (T) behaves as 



/2T27r2-3M2 

^(T) = J (34) 



Above Td, the pressure p = — / is given by 



15 4 167r2 



The internal energy is given by 



^ dp IGvr^T^ - 20M^T'^7r'^ - 15M^ 

£ = T—— — p = 7[ (36) 

dT ^ 807r2 ^ ^ 



The dimensionless interaction measure A, defined from the stress-energy tensor as A = 
i{TpT-{Ti:)o)/T\ is given by 

, e-3p M^iT^-Tj) 

Note for future reference that A falls off as l/T^ above T^. 

This first, analytically tractable example exhibits the one real shortcoming of these mod- 
els: the pressure can go negative at low temperatures, and shows non-monotonic behavior. 
A fully satisfactory theory would probably have the pressure identically zero in the confined 
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phase, or positive and very small if glueball effects were included. We could add a small 
constant to the free energy which would make p{Td) = 0, at the cost of slightly changing 
the graph of p/T"^ for temperatures just above T^. Such a term is similar to the appearance 
of the bag constant in the bag model, but with an important difference. In this model, any 
term independent of the temperature would also be independent of the Polyakov loop eigen- 
values, and would give the same contribution in both the confined and deconfined phases. 
The introduction of this new parameter, moreover, does not necessarily render the pressure 
physical in the region below T^. Since our intent is to model the behavior above T^, we view 
this behavior as a minor flaw. 



B. Model A for SU{3) 

In the case of SU{3), there are three eigenvalue — (p, and 0-1 = —0. At low 
temperatures, we have (p — 27i/3, giving Pp — 0. The free energy has the form 

fA = -T^^ + ^ [80^ (0 - + 0^ (0 - 2^f 

+ [2'/' (0 - tt) + (0 - 27r)] (38) 

As in the case of SU{2), a simple substitution is helpful. Deflning t/j — 27r/3 — 0, we may 
eliminate the linear terms and write /a in the form 

The presence of a ■0^ term indicates that the phase transition will be first order, as expected. 
The non-trivial minimum of /a, which represents the deconfined phase, is given by 



V'min = (40) 

The point at which ^l)■min develops an imaginary part is the spinodal point, given by = 
\/2M/t:. Below this temperature, the deconfined phase is no longer metastable. There 
is another spinodal point, associated with ■0 = 0. This temperature, which is given by 
Tl = 3M/2tt, is the temperature above which the confined phase is no longer metastable. 
A first order phase transition occurs when /(^/'min) = /(O), which gives the critical value 

^Jj = 2n/9 (41) 

and this in turn implies the deconfinement temperature is 

Td = —VlOM . 452 96M. (42) 

Comparing the two spinodal temperatures to Td, we have 

T 

-f . 993 81 

Td 

T' 

^ 1.0541 

Td 

which indicates a very narrow range of metastability around Td, consistent with a very weak 
first order phase transition. 
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V. MODEL B 



As will be demonstrated in section |VT|, model A gives a very reasonable approximation to 
the free energy and associated thermodynamic functions above Td, as obtained from lattice 
simulations. It is desirable to have a second model, to attempt to judge what features are 
universal. Our second model, like the first, is obtained by physically motivated considera- 
tions, but of a completely different type. 

Let us suppose that there is some natural scale R in position space over which color 
neutrality is enforced. In other words, net color is allowed in volumes of less than R'^, but 
the net color on larger scales is zero. We think of space as being divided up into cells of 
size R^, and assume each cell is large enough that the conventional density of states may be 
used. The general form of the partition function for a cell has been known for some time 
1T8I-E5I ; it has the form 



'.cell = exp 



PR' f cell] = J {d9) exp [-m'fpert {0)] (43) 



where {dO) denotes an integral over SU{N) Haar measure, and fpert{G) is as derived in 
section ||. Although f^eii turns out to have many of the properties associated with the 
deconfinement transition, it lacks a true phase transition, as any model based on integration 
over a finite number of variables must. In order to have a phase transition, there must be 
cell-cell correlations. We construct our model B by requiring that the eigenvalues 9 be the 
same in all cells. Such a condition might be derived by steepest descents, for example. This 
leads to 

fB {0) = Ucrt {0) - ^ In [J m (44) 

where J {0) is the Jacobian contained within the Haar measure {dO). For SU{N), its explicit 
form is 

up to a constant which is fixed by demanding that / {dO) = 1. At temperatures large 
compared to 1/R, minimization of Jb over the dj variables gives a free energy close to fceii- 

At low temperatures, the free energy is dominated by the measure term, which favors 
the same uniformly spaced pattern of eigenvalues found for model A. This behavior is very 
familiar from random matrix theory . At high temperatures, the free energy is dominated 
by fpert, and a phase transition occurs because of the conflict between the two terms. Phys- 
ically, the conflict is between minimizing the energy of the gauge boson gas and maximizing 
the entropy associated with color fluctuations. 

Model B clearly has features in common with model A. In both models, the correct 
perturbative behavior is built in, and the deconfining phase transition results from the 
interplay of two simple terms. Both models introduce a single mass scale which determines 
Td- However, the physical motivation for the two models is rather different, and the mass 
scale is introduced in different ways. 



11 



VI. SU{N) THERMODYNAMICS FOR N = 2,3,4, AND 5 



In this section, we present the pressure p, the energy density e and the interaction measure 
A as functions of T for SU{2) through SU{5) for both models. In the case of SU{3), we 
also compare the models with results from pure gauge simulations |^ . As discussed above, 
model A is analytically tractable in the cases of SU{2) and SU{3)] for the case of SU{4:) and 
higher, a numerical solution is easily obtainable. Model B, because it involves transcendental 
functions, must be solved numerically for all A^. 

For SU{2), both models predict a second-order phase transition, in agreement with 



universality arguments and lattice simulations [10|. For = 3, 4, and 5, the phase transition 



is first order in both models. Figures 1-3 show p, e, and A as a function of the dimensionless 
variable T/T^ for the case of SU{3). Results from lattice simulations, model A, and model 
B are shown. Since both models involve only a single dimensional parameter, these graphs 
have no free parameters. While neither model is in precise agreement with the simulation 
data, both are reasonable approximations over the range — 4.5 T^. The most notable 
discrepancy is that both models approach the blackbody limit faster than the simulation 
data. 

It is enlightening to plot p, e, and A, each divided by A^^ — 1, versus T/Td for = 2, 
3, 4, and 5. Figures 4-6 show the results for model A, and 7-9 for model B. Both models 
appear to be quickly approaching a finite large- A^ limit. 

Both the simulation data and the two models show power law behavior in A for suffi- 
ciently large T, consistent with A oc 1/T^. This is the asymptotic behavior found analyti- 
cally for SU (2) in model A. For the sake of a consistent analysis, we fit the simulation data 
and all models over the range ST^ — 4.5Trf to the simple power law behavior A oc l/T^. The 
results are presented in Table |. The results are compatible with 1/T^ behavior contaminated 
by subleading corrections. These results should be compared with the value p = 2.094(9) 
obtained by us from an analysis of the simulation data of [0. In fact, all of the curves can 
be fit extremely well over the range l.ST^ — 4.5Td by a curve of the form A = a/T"^ + b/T^. 

Considering that as formulated, models A and B have no free dimensionless parameters, 
they do a very good job of representing thermodynamic behavior. A more phenomenological 
approach would allow Mand R, respectively, to vary with temperature, allowing fits to 
lattice results. 



VII. QUARKS 

The effect of quarks can be straightforwardly included in an approximate fashion by 
adding the free energy of quarks propagating in a constant Polyakov loop background to the 
free energy of the gluons 0-^]. It is very important, however, to note that this procedure 
neglects chiral symmetry restoration. It is very likely that a unified model including order 
parameters for both deconfinement and chiral symmetry restoration is necessary to fully 



describe the quark-gluon plasma p8| , p9| . We limit ourselves here to a discussion of very 
heavy quarks and the leading order effect of light quarks. We will defer the subject of the 
interplay of deconfinement and chiral symmetry restoration to a later time. 

Quarks in the fundamental representation of the gauge group explicitly break Z{N) 
symmetry. As the quark mass m goes to infinity, this effect vanishes. On the basis of 
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theoretical models and lattice simulations, the expected effect of very heavy quarks is to 



lower the critical temperature |112| , p!3| . This line of first-order critical points in the T — m 
plane is expected to terminate in a second-order end point at some finite value of m. A 
low-temperature expansion |3^j3l|JT^ gives the free energy of a massive quark as 



U {0) = —r- E E —^K, {ndm) cos {nQ,) , (46) 

j=i „=i ^ 

a form suitable for (5m sufficiently large. 

We consider only the simple case of S'f/(3) using model A with a single heavy quark. 
Taking the parameter M = 1, the deconfinement temperature in the pure gauge theory 
is (m = oo) = 0.453M. The addition of a single heavy quark shifts the deconfinement 
temperature very slightly, with the first-order line terminating at a quark mass of m = 
2.57M. At this critical end point, the deconfinement temperature is Td {m = 2.57M) = 
0.974 Td (m = oo) = 0.441 M. 

For light quarks, the leading order behavior is proportional to and independent of 
the quark mass, and therefore independent of chiral behavior. It is given by [^-Q 

f (^) = Nf (-^ - —tfj - -i? + _ 1 y4 (47) 

^\ 4860 81 ^ 9^ 97r^ 67r2^; ^ ' 

where N j denotes the number of light flavors. Recall that the conflned phase of the pure 
gauge theory was given by ■?/' = 0. The term linear in appearing in (■?/;) directly indicates 
the breaking of the 2'(3) invariance of the pure gauge theory. Detailed analysis shows that 
with N f = 2, the deconflnement transition is replaced by a smooth crossover. The only scale 
in this model is again the parameter M, which need not have the same value found in the 
pure gauge theory, and a rapid rise in all thermodynamic quantities begins at a temperature 
around M/4. 



VIII. CONCLUSIONS 

We have developed two phenomenological equations of state for the quark-gluon plasma, 
which reproduce much of the thermodynamic behavior seen in lattice simulations. These 
models have many attractive features. Not the least of these is simplicity. Both models 
introduce a single new parameter. For the most important case of SU{3), the free energy 
is obtained by minimizing over a single variable. Both models correctly predict the order of 
the deconflning phase transitions for SU{2) (second order) and SU{3) (weakly flrst order). 
Both models predict flrst order transitions in SU{4) and SU{5) as well. The pressure and 
other thermodynamic quantities vary rapidly in the range Td — 5Td in both models, and both 
appear to have smooth large N limits. 

Although the two models have very different phenomenological origins, the numerical 
value of the parameters introduced are reasonable for the case of SU{3). If we take the 
deconflnement temperature in a pure SU{3) gauge theory to be 270 MeV, then model A 
gives a value for M of 596 MeV. This is a plausible value for a constituent gluon mass. In 
model B, we flnd that i?is 1 fermi, which is of course a typical hadronic scale. Thus in both 
models the phenomenological parameter introduced has a reasonable value. 
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It is useful to compare these models to the naive Bag model equation of state, which is in 
common usage in phenomenological applications [Q. The Bag model also introduces a single 
dimensional parameter, the bag energy density B, but gives a first order transition for all 
N. The Bag model pressure approaches the high-temperature, blackbody limit faster than 
both model A and model B, and also lattice results. The graph of e/T^ versus Tfor the Bag 
model is monotonically decreasing, in complete disagreement with lattice results. Finally, 
the Bag model predicts a behavior for A, which is ruled out by lattice results. As we 
have seen, model A and model B both show a behavior in A, which is compatible with 
SU{3) lattice data. 

As we have formulated them, both model A and model B have no free parameters 
once Td is fixed. It is clear that by allowing the parameters M and R to depend on the 
temperature, a better fit to lattice data can be obtained at the cost of introducing additional 
phenomenological parameters. We plan to explore this issue for the most important case of 
SU{3). 

The success of these phenomenological models strongly suggests a theoretical point of 
view on the nature of the deconfinement mechanism, independent of the specific gauge theory 
under study. In confining theories, such as pure gauge theories, the eigenvalue distributions 
of the Polyakov loop are peaked at low temperature around values evenly spaced about the 
unit circle in such a way that the expectation value of the Polyakov loop is zero. Such 
behavior is characteristic of random SU (N) matrices. At the deconfining transition temper- 
ature, the peaks of the eigenvalue distributions change, and move towards 1, which is the 
asymptotic limit as T goes to infinity. The consistency of this picture can be tested using 
lattice Polyakov loop eigenvalue distributions for SU{4) and higher. In theories with light 
quarks, the deconfinement transition may be replaced by a rapid crossover, which is again 
associated with the motion of the Polyakov loop eigenvalues. We believe that this picture of 
deconfinement coupled with a field-theory inspired model of chiral symmetry breaking has 
the potential to fully model the equation of state of the quark-gluon plasma. 

Associating deconfinement with changes in Polyakov loop eigenvalue distributions gives 
some perspective on recent attempts to model the plasma equation of state using the hard 
thermal loop (HTL) approximation p^-|35|. The HTL approximation is a resummation 
of high-temperature perturbation theory, and may provide a good explanation of the slow 
asymptotic approach to the blackbody limit at very high temperatures. However, if the 
picture we have advocated here is correct, then no approach based on perturbation theory 
can explain the behavior of the plasma in the crucial region from to ST^, because the 
Polyakov loop eigenvalues are assumed to be 1. In fact, preliminary numerical simulations 
of pure SU{3) gauge theories using Nt = 4 lattices indicate that the eigenvalues of indi- 
vidual Polyakov loops move much more slowly towards 1 than either model predicts. It is 
possible that a hybrid approach combining hard thermal loops with the minimal amount of 
phenomenology used in our models would be quite successful in reproducing the quark-gluon 
plasma thermodynamics over the entire range of temperatures above T^. 

Finally, there is the question of the connection of these models to candidate explanations 
for confinement having a deeper basis in the underlying gauge theory dynamics. Both model 
A and model B attribute confinement to a single, simple term in the free energy. Derivation 
of a similar term from one of these more fundamental candidate explanations of confinement 
would be very satisfying. A proper field-theoretic basis for the motion of the Polyakov loop 
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eigenvalues would answer many questions. 
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TABLES 



TABLE L The parameter p obtained from the fit A cx 



Model A N=2 L833(3) 

Model A N=3 L878(2) 

Model A N=4 L889(2) 

Model A N=5 L893(2) 

Model B N=2 2.377(4) 

Model B N=3 2.398(3) 

Model B N=4 2.436(3) 

Model B N=5 2.457(3) 
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